
' Figure 1. (a): A cross-cut along the mid-plane of the disk, z = 0, at t = 10.4 Myr. Density contour lines and magnetic field lines (running 
horizontally) are plotted. The numbers near the bottom-right corner represent, respectively, the time of the snapshot, the maximum and 
minimum levels of the contours, e.g., 21 contour lines are plotted for the values from 10 _310 cm — 3 to 10 — 38 cm -3 , (b): The same as (a) 
but for the x = plane. Density contour lines and velocity fields are plotted. The magnetic field is running perpendicular to this plane 
and its distribution is represented by dotted contour lines. The maximum speed plotted in the panel is also shown near the bottom-right 
corner, (c): The same as (b) but for the y = plane. 
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Figure 2. Physical quantities are plotted along the x— (a), y— (b), and z— axis (c). The density n(cm -3 ) (solid lines), the temperature 
T(10 4 K) (dotted lines), the magnetic flux density B x ([iG) (dashed lines), and the velocities along the respective axes |t)|(km s _1 ) 
(long-dashed lines) are plotted. This is a snapshot at the age of t = 10.4 Myr. 
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Figure 3. The same as Fig.l, but for the age of t = 22 Myr. 
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Figure 5. Expansion law of a hot cavity formed in a superbubble. The size of the hot cavity is measured with (x c ,y c , z c ). Models A and 
B correspond to the cases that a superbubble is formed in ISM with a uniform magnetic field distribution. In Models C, D, E, and F, 
the magnetic field strength was assumed to decrease upwardly, while the Alfven speed S/(47rp) 1 / 2 is constant. Model Q, non-magnetic 
superbubble, is also plotted. 
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Figure 6. A supcrbubble driven by an active OB association with a luminosity 10 times larger than the previous model. This is a 
snapshot at the age of t = 20 Myr. 
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Figure 7. The same as Fig.l but for the model of Bo oc p 1 / 2 . Panels (a), (b) and (c) are the snapshots of t = 10 Myr. Panels (d), (c) 
and (f) are those of t = 22.8 Myr and (g), (h) and (i) arc for t = 30 Myr. 
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Figure 7. continued. 
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Figure 7. continued. 
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Figure 8. The same as Fig.7 but for a model with a less mechanical luminosity of Lsn = 10 37 erg/s. The snapshot corresponds to t = 40 
Myr. 
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Figure 9. Expansion law for a thin-shell model. External interstellar pressures are chosen equal to p ut = 1.7 X 10 — 12 erg cm -3 (solid 
line), pout = 1 X 10 -12 erg cm~ 3 (dotted line), and p ut = (dashed line). The Weaver et al.'s (1977) solution is similar to a dashed 
curve and the point when the ram pressure is equal to the interstellar pressure with p ou t = 1.7 X 10~ 12 crg cm~ 3 is also plotted by a 
asterisk. The equivalent radii which are defined as r equ i v = (3V ho t/47r) 1 / 3 for Models A and C are also illustrated with dash-dotted lines 
(lower:A; uppcr:C). 
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ABSTRACT 

The importance of the interstellar magnetic field is studied in relation to 



' the evolutions of superbubbles with a three-dimensional (3D) numerical mag- 

^vq ' netohydrodynamical (MHD) simulation. A superbubble is a large supernova 

o ' 
o 

oo 

G\ ■ the size reaches 2-3 times the density scale-height, the superbubble expands 



remnant driven by sequential supernova explosions in an OB association. Its 
evolution is affected by the density stratification in the galactic disk. After 



preferentially in the z-direction. Finally it can punch out the gas disk (blow- 
out). On the other hand, the magnetic field running parallel to the galactic 
disk has an effect to prevent from expanding in the direction perpendicular to 



CO 

the field. The density stratification and the magnetic fields have completely 



opposite effects on the evolution of the superbubble. We present results of 
3D MHD simulation in which both effects are included. As a result, it is con- 
cluded that when the magnetic field has a much larger scale-height than the 
density, even for a model that the bubble would blow out from the disk if the 
magnetic field were absent, the magnetic field with the strength of 5 /iG can 
confine the bubble in \z\ ^ 300 pc for ~20 Myr (confinement). In a model 
that the field strength decreases in the halo in proportion to B oc p 1//2 , the 
superbubble eventually blows out like a model of B = even if the magnetic 
field in the mid-plane is as strong as B = 5/iG. 

Key words: superbubbles - supernova remnant - magnetic fields. 
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1 INTRODUCTION 

A superbubble is a complex consisting of an OB association, surrounding X-ray emitting 
hot gas, and a corresponding HI hole/shell. Three examples are well-known in our Galaxy: 
Cygnus (Cash et al. 1980), Orion- Eridanus (Cowie, Songaila & York 1979; Reynolds & Ogden 
1979), and Gum nebula (Reynolds 1976). Superbubbles are found as HI shells and holes in 
external galaxies, such as LMC (Meaburn 1980; Dopita, Mathewson & Ford 1985), M31 
(Brinks & Bajaja 1986), M33 (Deul & den Hartog 1990), M101 (Kamphuis, Sancisi & van 
der Hulst 1991), and so on. Since the sizes of these objects are in the range of 100 pc - lkpc, 
this can not be explained by a single supernova explosion [the size of an ordinary supernova 
remnant (SNR) is ^ 50pc]. The amount of energy required for such a superbubble reaches 

5 x 10 51 erg — 10 54 erg ( Tenor io-Tagle & Bodenheimer 1988). There are two models proposed 
for formation of superbubbles: (1) a large SNR driven by sequential supernova explosions in 
an OB association and (2) a complex formed by a collision of a high-velocity cloud and the 
galactic disk (Tenor io-Tagle 1991). Here, we confine ourselves to the first model and discuss 
the evolutions. For review papers of this field, see Tenorio-Tagle & Bohdenheimer (1988), 
Spitzer (1990), Tomisaka (1991), Bisnovatyi-Kogan & Silich (1995). 

Here, we summarize the evolution very briefly. After the size of bubble exceeds the density 
scale- height, the bubble becomes elongated in the direction perpendicular to the galactic 
disk (Tomisaka & Ikeuchi 1986; Tenorio-Tagle, Bodenheimer & Rozyczka, 1987; MacLow 

6 McCray 1987). It is shown that when the mechanical luminosity released by sequential 
supernova explosions is high, the expansion of the bubble to the halo is accelerated and the 
hot gas contained in it flows into the galactic halo (galactic fountain), finally. However, the 
magnetic fields running parallel to the galactic disk prevents the gas from flowing in the 
vertical direction (perpendicular to the fields). 

Tomisaka (1990) studied the adiabatic evolution of a superbubble in uniform magnetic 
fields and showed that the bubble driven by a mechanical luminosity of Lsn — 3 x 10 37 ergs _1 
is confined in the galactic disk by the effect of the magnetic field provided its strength is as 
large as B = 5/zG. This shows the superbubble is confined in the galactic disk, if (1) the 
magnetic fields have a large scale-height, H B , and (2) L SN ~ 3 x 10 37 ergs _1 and B <; 5/iG. 
However, this result may be affected by the assumption of the adiabatic gas. The interstellar 
magnetic fields seems to play an important role in the evolution of a superbubble (see also 
Mineshige, Shibata, & Shapiro 1993, Ferriere, MacLow, & Zweibel 1991). 
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In the present paper, a full magnetohydrodynamical calculation has been done including 
the radiative cooling. Further, the effect of distributions of magnetic field strength is studied. 
Plan of the present paper is as follows: in section 2 is given an analytical estimate to a 
threshold mechanical luminosity under which the superbubble is confined in the galactic disk. 
From this estimation, we choose the values for the parameters Lsn ; Bq, and Hb- Numerical 
method and other assumptions are also presented in §2. Section 3 is for the numerical 
result, in which the effect of the magnetic fields is shown. In section 4, we will discuss the 
observability of the superbubble in external galaxies. Using the evolution obtained, we will 
show that a large fraction of the interstellar space is occupied with superbubbles. 

2 ANALYTICAL MODELS 

As a first step, the expansion law of a superbubble is studied with a simple approximation of 
spherical symmetry. [Historically, it was done by Bruhweiler et al. (1980) and Tomisaka, Habe 
& Ikeuchi (1981).] A self-similar solution for a wind-blown bubble by Weaver et al.(1977) 
shows us the expansion law for the shock wave driven by a steady energy release as 

R s = 271pc(L SN /3 x 10 37 erg s" 1 ) 175 

x (p /6 x l(T 25 g cm- 3 )- 1/5 (t/10Myr) 3/5 , (1) 

where Lsn is a mean mechanical luminosity ejected from supernovae, i.e., Lsn = Eq/At 
using the total energy ejected from an SN (Eq) and the mean interval of two SNe (At). 

Since a typical size of the superbubble (R s ) is much larger than a density scale-height 
of the galactic disk, H = / °° p(z)dzj p(0), the evolution should be affected by the density 
stratification in the z— direction. This problem is now two-dimensional. A frontier of this 
problem was Chevalier & Gardner (1974), who studied the evolution of an SNR in the galactic 
halo using a 2D hydro-code, although a small number of zones could be used at that time 
as 20 x 39. This kind of problem became attacked seriously after mid-1980's using large 
2D hydro-codes (Tomisaka & Ikeuchi 1986, 1988; Tenorio-Tagle, Bohdenheimer & Rozyczka 
1987, Tenorio-Tagle, Rozyczka & Bohdenheimer 1990; MacLow, McCray & Norman 1989; 
Igumentshchev, Shustov & Tutukov, 1990) and semi-analytic approximations (MacLow & 
McCray 1987; Bisnovatyi-Kogan, Blinnikov & Silich 1989; Koo & McKee 1990; 1992). 

Summarizing results of the numerical studies above, it is concluded that in the first 
stage, the expansion is spherical as long as the size is much smaller than the density scale- 
height. However, when the Lsn is sufficiently large, a bubble enters a new phase, in which 

© 1997 RAS, MNRAS 000, §-0 



4 K. Tomisaka 

an expansion in the direction of the density gradient becomes accelerated. In an exponential 
atmosphere as p = po exp(— \z\/H), Mac Low et al. (1989) showed that the shock front 
is accelerated after it passes z ~ 2.9H. Although the exact numerical factor is dependent 
on how the density distributes, p(z), the numerical simulations insure that the shock is 
accelerated upwardly, after the size is larger than (2 — 3) x H. Finally, the superbubble 
punches out the gas disk and a hot gas contained in the bubble flows into the halo of the 
galaxy. 

How about a less-energetic bubble? If the thermal pressure po of the ambient interstellar 
medium (ISM) is dominant over the post-shock ram pressure, i.e., pov 2 ^ po, before the 
break-out, the bubble seems to be confined in the disk. Using equation (1), the ram pressure 
becomes equal to the thermal pressure at the age of 

t P ~ 79MyrL^ 2 ^ /4 P :^ 4 , (2) 

where L 38 = L SN /10 38 erg s _1 , p_i 2 = £>o/10~ 12 dyn cm -2 , and n represents the number 
density of ISM. We assume here that the condition for the break-out is R s (tp) > aH (Koo 
& McKee 1992). Equating R s (tp) = aH, we obtain a threshold luminosity above which the 
break-out occurs as: 

L crit = 0.59 x 10 37 erg sr\o? /h)I%nZ V VS, (3) 

with H 2 = if/lOOpc. The value of a ~ y/E is taken from an estimate by numerical results by 
Mac Low & McCray (1988). Thus, if is much smaller than L crit , the bubble is eventually 
confined in the galactic gas disk. 

2.1 Effect of the interstellar magnetic fields 

The magnetic pressure in ISM is expected as 

Pmag = B 2 /8n ~ 10- 12 dyn cm- 2 (B /5fiG) 2 . (4) 

Taking both magnetic (p ma , g ~ 10~ 12 dyn cm -2 ) and thermal pressures (pth = nkT ~ 
2400A; ~ 7 x 10 _13 dyn cm -2 ) into account in equation (|^), the critical luminosity becomes 

L crit ~ 3 x 10 37 erg s' 1 (iJ/180pc) 2 (no/0.3cm- 3 )- 1/2 

x (po/1.7 x 10" 12 dyn cm^ 2 ) 3/2 , (5) 

where we take a typical value for H as 180 pc according to a model by Dickey & Lockman 
(1990). This shows that the magnetic field seems to play a role for a superbubble with Lsn — 
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3x 10 37 erg s _1 . This corresponds to a mean interval between subsequent supernova explosions 
in an OB association of At ~ 10 6 yr and this is not inconsistent with that estimated from SN 
rate in our galaxy (Tomisaka et al. 1981). That is, the superbubble might be magnetically 
pressure-confined in the disk. 

Another possibility for a superbubble to be confined in the disk is that the mechan- 
ical luminosity decreases before the shock front is accelerated in the z— direction. Since 
the age when the size R s exceeds aH is equal to r acc ~ 6MjT(a/y/5) 5 / 3 H^ 3 (Ls^s/3 x 
10 37 erg s _1 ) _1 / 3 (p/6 x 10~ 25 g cm" 3 ) 1 / 3 , in the case that the lifetime of the active phase 
of an OB association r\a e is shorter than r acc , the expansion is stalled due to the decrease 
of the mechanical luminosity. If we assume the sequential supernova explosions continue at 
least riifc <: 40Myr (e.g. Shull & Saken 1995), only for an extremely poor OB association as 
L SN < 1 x 10 35 erg s _1 (ri ifc /40Myr)~ 3 (a/v / 5) 5 -ff|(p/6 x 10~ 25 g cm" 3 ) the superbubble burns 
out before the shock is accelerated upwardly in the z— direction. 



2.2 Models and basic equations 

We choose parameters for interstellar medium as follows: assuming the disk-plane (z = 0) 
number density no = 0.3cm -3 , the gas temperature To = 8000K, these give a pressure on 
the disk-plane as po ~ 7 x 10~ 13 dyn cm" 2 . The density distribution perpendicular to the 
disk is taken based on a model proposed by Dickey & Lockman (1990) as 

, 2" 



n(z) 



0.395 exp 



0.566 

+0.107 exp 
+0.064 exp 



1 

2 V90pc 



225pc 
1*1 



403pc 



(6) 



this yields an effective scale-height as H e{i = J °° pdz/p ~ 180pc. Applying equation (^j), 
model parameters chosen here give a critical luminosity of L cr - lt ~ 3 x 10 37 erg s _1 . Parameters 
taken are summarized in table 1. As an initial condition, we began with a hydrostatic state 
such as the thermal and magnetic pressures balance with the gravity due to the galactic disk 

as 

1 d(p th + B 2 /8n) 

9z{z) = j. . (7) 

p oz 

Actually, we assumed the distributions of gas and magnetic pressures first. Then, we chose 
the gravity consistently. In addition, we add a halo with constant density for z > 0.8 kpc 
to avoid numerical difficulties. The gravitational acceleration assumed here is somewhat 
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Table 1. Model Parameters 



Model 


no -Bo 


T 


L SN B x {z) ifinalF] 


-f^cquiv[j 




(cm" 3 ) (pG) 


(K) 


(erg s" 1 ) (Myr) 


(PC) 



A 0.3 5 8000 3 x 10 37 const 22 255 

B 0.3 5 8000 3 x 10 3S const 20 567 

C 0.3 5 8000 3 x 10 37 oc p 1 / 2 36 592 

D 0.3 3 8000 3 x 10 37 oc p 1 / 2 30 563 

E 0.3 5 8000 1 x 10 37 oc p 1 / 2 40 196 

F 0.3 5 8000 3 x 10 38 oc p 1 / 2 11 470 

Q 0.3 8000 3 x 10 37 no B 50 1150 



a Numerical run continues up to this age. 
6 Equivalent radius, eq. (|l9|), at tfl Ila i. 

different from that is derived from observations of distributions and motions of K giants 
(Spitzer 1978). For example, equation (7) gives 7 x lCT 9 cm/s 2 for z = 200 pc and 4 x 
10 _9 cm/s 2 for z = 400 pc, while the values derived from K giants for respective heights 
are equal to 4 x 10 _9 cm/s 2 and 6 x 10 _9 cm/s 2 . Since the gravitational free-fall time scale 
~ (z/g z ) l l 2 ~ 20Myr(z/lkpc) 1//2 (g z /8 x 10~ 9 cm/s 2 ) -1 / 2 is comparable to a typical age of 
the superbubble, the gravitational acceleration may play a role in the late phase and then 
the gravitational acceleration assumed here may underestimate the effect of gravitational 
confinement of the superbubble. 

The magnetic scale-height of our Galaxy can be estimated with the rotation measures of 
extragalactic radio sources (Spitzer 1978). This gives a scale- height of a thicker component 
of the magnetic disk whose e-folding scale-height is estimated as 1.2 ±0.4 kpc (Han & Qiao 
1994 and references therein). First, in Models A and B we assume the strength of magnetic 
fields is constant (B (z) =const) irrespective of the height from the disk midplane, in which 
we will see the effect of this thick magnetic disk. In addition, from the rotation measure 
observations of galactic pulsars, it is shown there is an additional narrower component of 
the magnetic disk whose Gaussian scale- height is as large as 170 pc (Thompson & Nelson 
1980). To take this narrower component into consideration, we assume the magnetic field 
strength decreases in proportional to the square root of the gas density in Models C-F. 

Basic equations to be solved are the ideal MHD equations as follows: 

f + V-pv = 0, (8) 

dpv 1 

+ V • pvv = -Vp + pg + —V x B x B, (9) 
ut 4tt 

g + V • ev = -A - pVv, (10) 
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with 

7-1 

f)~D 

= V x (v x B), (12) 

where A(T, p) = L(T)n 2 represents the radiative cooling. Here, we fitted values derived by 
Raymond, Cox & Smith(1976) with a third-order spline function. Other symbols have usual 
meanings. Equations (||), ©, ( |10|) and ( |12|) represent, respectively, the continuity equation, 
the equation of motion, the energy conservation equation, and the induction equation under 
the ideal MHD condition (with infinite conductivity). The adiabatic exponent 7 is taken 
equal to 5/3. 

Numerical method we employ is a finite-difference three-dimensional magnetohydrody- 
namics code. Specifically, a "monotonic scheme" (van Leer 1977; Norman & Winkler 1986) 
is adopted to solve the hydrodynamics equation, equations (|8|), (|) and (pi]). On the other 
hand, a "constrained transport" method (Evans & Hawley 1988) and a "method of char- 
acteristics" (Stone & Norman 1992) are used to solve the magnetic part of equation @ 
and the induction equation (|12|). The code is well vectorized and parallelized. Thus, it is 
executed very efficiently on a vector-parallel machine Fujitsu VPP300 with 16 processors, 
e.g., it needs only 0.9 sec - 3.8 sec to proceed one time step for grids 161x161x261 - 
271x271x451, respectively. The number of zones is increased just before the superbubble 
expands and it touches one of the numerical boundaries. The grid spacing (Ax, Ay, Az) is 
chosen spatially constant as 5pc. A superbubble is located at the origin x = y = z = and 
energy and mass are ejected in a number of cells near the origin every time step At to the 
amount of Lg^At and MsnA£. The mass ejection rate Msn is chosen as 3 x 10 22 g s _1 for 
-^sn = 3 x 10 37 erg s _1 . A superbubble is assumed symmetrical with respect to x = 0, y = 
and z = planes. Thus, 1/8 of the volume of the bubble is actually calculated. Correspond- 
ing three symmetric boundary conditions are set for x = 0, y = 0, and z = planes. The 
other outer boundary conditions are set as free boundaries. 



3 NUMERICAL RESULTS 

3.1 Model with uniform magnetic fields 

First we compare two models with Lsn < -^crit and Lsn > -^crit- Next, the effect of the 
distribution of magnetic field strength is studied. In Model A, uniform magnetic fields are 
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assumed. Figure 1 shows the cross-cuts of the superbubble at the age of 10.4 Myr after the 
first SN exploded. A non-magnetic superbubble consists of two parts, that is, an inner cavity 
which is occupied with a hot gas ejected by SN explosions and a shell which consists of ISM 
compressed by the effect of inner high-pressured cavity. This basic structure is common in 
magnetized superbubbles. Figure 1 (a) indicates that the bubble elongates to the direction of 
global magnetic fields on the mid-plane of the galactic disk (z = plane). Hot low-density 
gas ejected from SNe fills a hot inner cavity which extends up to x c ^ 220 pc (near the 
x— axis) and y c ^ 150 pc (near the y— axis). Outside the hot cavity an outward-facing fast- 
magnetosonic shock front is propagating, which is traced as a series of sharp bends of the 
magnetic fields. Since the magnetic fields bend as they reach the front after passing the 
front, this is a fast-mode shock front. The structure perpendicular to the disk is seen in Figs 
1(b) and 1(c). In Fig. 1(b), iso-density lines (solid lines) and contour lines of magnetic flux 
density (B x ) in the direction perpendicular to the plot plane (e.g., x = plane) are shown. 
The outer MHD shock front is clearly traced in these figures both by the distribution of the 
magnetic flux density (dotted lines in Fig. 1(b)) and bends of the iso-density contour lines 
(solid lines in Figs 1(b) and (c)). It should be noticed that the shock front is detached from 
the hot cavity and it runs much faster than the contact surface between the ejected matter 
and accumulated ISM. This indicates that a thick shell is formed between the contact surface 
and the fast-mode shock front. This is much different from the non- magnetic superbubble 
(Tomisaka & Ikeuchi 1986; MacLow et al. 1989), in which a hot cavity is surrounded by a 
cooled thin shell. This is understood as follows: since the magnetic force directs perpendicular 
to the fields, compression perpendicular to the field is blocked and the shell becomes much 
thicker than that of a non-magnetic superbubble. 

In Fig. 2, physical quantities are plotted along the axes. Temperature distribution (dotted 
lines) clearly indicates the boundary between a hot cavity with T ~ 2 x 10 6 K and a shell 
with T ^ 10 4 K. Since the hot cavity is filled with a matter ejected by SN explosions, in the 
cavity there is no magnetic field which has its origin in the interstellar space. It is shown 
that on the thin shell is formed around x ~ 250 pc. A magnetosonic wave front 

is propagating at y s ~ 350 pc on the y— axis and at z s ~ 600 pc on the z— axis, while the 
size of the hot cavity is no more than y c ~ 150 pc and z c ~ 200 pc. In a region y c < y < y s 
and z c < z < z s , a thick shell is expanding. After 15 Myr, another magnetosonic wave is 
detached from the the contact surface and propagates upwardly inside the thick shell. 

In t = 22 Myr, the hot cavity expands up to x c ^ 420 pc, y c ^ 160 pc, and z c ^ 330 pc 
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(Figs 3 and 4). The elongated structure to the x— direction is characteristic of the magnetized 
superb ubble. The magnetosonic wave front seen in Fig. 1(b) has propagated away from the 
frame of Fig.3(b). Contraction to the contact surface from outside occurs near the mid-plane 
of the disk [v x is negative in the region of 450pc ^ x ^ 650pc in Fig. 4 (a) and v y is negative 
in the region of 150pc ^ y ^ 250pc in Fig.4(b)]. Since the temperature of this region is lower 
than the average of the shell, this contraction seems to be driven by an effective cooling near 
the contact surface. 

Although the fast-mode magnetosonic wave has propagated away out of the panel of Fig. 3 
(z > 1.3 kpc), a hot cavity is confined below z ^ 320 pc. Expansion of contact surface in each 
direction, x c , y c , and z c , is plotted against the age in Fig. 5. Compared with Model Q (non- 
magnetic case), it is shown that in the x— direction the hot cavity expands faster than Model 
Q. This means that the pressure in the cavity which drives the outer shell works efficiently 
near the x— direction, since there is no magnetic force in this direction. In contrast, in the 
y— direction the expansion is suppressed due to the magnetic tension force. The expansion in 
the y— direction ends before t ^ 17Myr, while the expansion parallel to the global magnetic 
fields continues as long as the simulation continues for t — 22 Myr. Contraction in the 
y— direction seems to be driven mainly by the Lorentz force which works to compress the 
hot cavity in the negative y— direction. The expansion speed in z— direction is still positive 
at the age of t = 22Myr but it is heavily decelerated by the effect of magnetic tension force 
working downward. The height of the non-magnetic superbubble in Fig. 5 reaches z c ~ 770 
pc in t = 22Myr and it indicates only a small deceleration. While, a superbubble with the 
same energy release rate and age but formed in a magnetized ISM with B = 5/xG is almost 
confined to the galactic disk z c ~ 300 pc. From this model, it is shown that the hot gas 
contained in the cavity does not break through the disk and is confined in the disk. Although 
the magnetic force works downwardly, we found no down-flow in the thick shell. 

3.2 The effect of the supernova rates 

In this section, we compare the models with different supernova rates. That assumed in 
Model B (Lsn = 3 x 10 38 erg s -1 ) is 10 times larger than that of Model A. In Fig. 6 the 
snapshot at t — 22Myr is plotted. Compared with Figs 3 and 4 (t — 22Myr), it is clear 
that the superbubble expands much faster than Model A. The top of the hot cavity reaches 
z c ~ 850pc, while in Model A it is only 330 pc. Figure 6(a) indicates that cross-cut view 
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along the surface of z = is rounder than that of Fig. 3 (a). This seems to come from a fact 
that the magnetic force becomes relatively unimportant compared with the ram pressure 
driven by the sequential explosions in this case. Comparing Figs 6(b) (x — plane) and 
6(c) (y = plane), it is shown that the bubble is expanding preferentially in the direction 
parallel to the magnetic fields. Since the bubble is confined by the magnetic fields in y— 
and z— directions, the bubble expands along the magnetic field lines after reaching the halo 
(400pc £ z £ 600pc). 

Figure 6(b) clearly indicates the bubble consists of two blobs which are connected ver- 
tically with each other. This is also seen in other models. Although the bubble is almost 
spherical initially, the expansion of the top of the shell is accelerated progressively due to a 
steep vertical density gradient and the shape is deformed. 

Shown in Fig.5(c), the expansion in the z— direction, z c , is decelerated due to the mag- 
netic force but dz c /dt is still positive. Although the expansion speed in the y— direction 
becomes negative (contracting) in a later phase of Model A, this model shows a positive 
velocity (expanding) throughout the evolution. This is also explained by the fact that the 
ram pressure force is relatively important against the magnetic tension which works to de- 
celerate the expansion in the y— direction. In the x— direction, expansion continues without 
any deceleration. 

Summarizing the results for Models A and B, the expansion in the z— direction is more 
or less suppressed by the uniformly distributed magnetic fields running in the a;— direction. 
The hot gas seems to be confined by a uniform magnetic field with the strength of Bq = 5//G 
to the extent of \z\ ^ 400pc (Model A) and \z\ & lkpc (Model B), when the scale-height of 
the distribution of magnetic field strength is much larger than that of the density as Models 
A and B. These models correspond to the infinite Hb- In the following section, we focus on 
the case with a finite Hb- 

3.3 The effect of distribution of the magnetic fields 

In this section, we focus on the effect of distribution of the magnetic fields. In Model C, 
we assume that the strength of magnetic fields decreases from disk to halo. The ratio of 
the thermal energy density to the magnetic one is assumed constant. That is, since the gas 
temperature is constant, the magnetic field strength changes according to the density as 




(13) 
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where B and po represent, respectively, the magnetic field strength and the gas density 
on the mid-plane of the galactic disk (z = 0). This gives a scale- height measured by the 
magnetic energy 



equal to the density-scale height H. The density distribution is again assumed as a model 
proposed by Dickey and Lockman (1990). Then, since the gas is supported by the magnetic 
pressure in contrast to Models A and B, the corresponding gravity becomes weak. Model C 
assumes the same mechanical luminosity Lsn as Model A. 

The snapshots at the ages of t = lOMyr (panels a-c), t = 22.8Myr (panels d-f), and 
t = 30Myr (panels g-i) are plotted in Fig. 7. As shown in Fig. 7(b), the hot cavity consists 
of two blobs similar to Figs 3(b) and 6. This panel captures a structure just after a top of 
the spherical shell is deformed and hot gas outflowing through the hole is re-expanding into 
the halo region. The reason why the hole is made seems to an instability which destructs 
the shell. 

The pressure in a hot cavity decreases with time as from p ~ 6000 K cm -3 at the age 
of t — lOMyr to p ~ 700 K cm -3 at t — 30Myr. The volume of the hot cavity increases 
more than 20 times in this time span. If the bubble were to expand exactly in an adiabatic 
fashion, the pressure would decreases much, ~ 1/150. However, it is certain that a rapid 
expansion into the low-density halo reduces the pressure rapidly after t = 10 Myr. 

In this model, as increasing the height, the effect of magnetic fields decreases. Thus, after 
the bubble escapes from the disk, it becomes hard to confine the hot gas with the magnetic 
fields. Comparing with Fig. 3, even with the same mechanical luminosity as 3 x 10 37 erg s _1 , 
it is evident that the superbubble is expanding faster than that of Model A. 

Comparing Fig.7(d) with Fig.3(a), as for the mid-plane (z = 0) the area occupied with 
a hot gas is smaller than that of Model A. Contraction driven by a magnetic tension occurs 
also in this model. In Fig.5(b), the time evolution of y c is illustrated. This shows that the 
maximum size of the superbubble on the mid-plane perpendicular to the magnetic field, 
y c , decreases after t ~ lOMyr (Figs 7(a)-(c)). This occurs after the pressure of the cavity 
decreases due to a rapid expansion of the bubble into the halo. The shell moves reversely 
in the y— direction by the effect of the magnetic tension. In the course of the contraction, a 
part of the shell in y > and that in y < merge together. At that time, the size in the 
direction parallel to the magnetic field, x c , decreases abruptly (t ~ 22Myr). 
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Figures 7(b), 7(e), and 7(h) indicate that the hot cavity is surrounded by a thick shell 
viewing from the direction of the magnetic field. At the age of t = 30 Myr, in a region 
near the z— axis (z ~ 200 — 300 pc), the hot matter ejected with supernova explosions flows 
through a narrow channel like a throat. A long sheet is extending from the throat, which 
was originally a shell on the top of the bubble. Figures 7(c), 7(f), and 7(i) shows that the 
shell propagating along the magnetic fields is also thick. The interface between the shell and 
the hot interior is clearly seen as a sharp density discontinuity. 

It is concluded that in a galactic disk whose magnetic scale-height, Hb, is as small as 
that of the density, H, a superbubble evolves like a bubble in a water. That is, a low density 
bubble rises by the effect of buoyancy. 

3.4 Low luminosity bubble 

In the preceding subsection, it is shown that an outflow of hot gas is realized in a galactic 
disk with a small magnetic scale-height. Then, does a less luminous bubble make the hot gas 
flow into the halo even from the same initial conditions as Model C? In Model E, a poor OB 
association is assumed with a mechanical luminosity of Lsn = 10 37 ergs _1 . The expansion 
law for this model is also plotted in Fig. 5. Owing to a low luminosity, the expansion is slower 
than Model C. 

Figure 8 shows the snapshot at the age of t = 40Myr. As shown in Fig. 5, the snapshot 
corresponds to the structure after the contraction occurs due to the magnetic tension. Com- 
paring with Model C, it is evident that the volume occupied with a hot cavity is small. A 
fraction of the hot cavity to the entire volume inside the magnetosonic shock wave front is 
also much smaller than that of Model C. From this figure, since the volume of hot gas is small 
and the height of the top z c is no more than ~ 500 pc, it is concluded that less luminous 
bubbles do not play important roles as an origin of the interstellar hot gas of T ^ 10 6 K. 

4 DISCUSSION 

4.1 Blow-out or confinement? 

As shown in the preceding section, although in Model B the mechanical luminosity is much 
larger than the critical luminosity of equation^), the vertical expansion is much decelerated 
by the effect of magnetic fields. This shows that L SN < L crit may be only a sufficient condition 
for confinement of the superbubble. 
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The expansion of a spherical shell driven by a pressure in the hot cavity p is formulated 
as follows: 

dMv 9 . , , 

— = 4nR 2 (p-p out ), (15) 

— = 4nR 2 vp out (16) 

^ = L SN -A7rR 2 pv, (18) 

where, M, i?, t>, p out , p ou t and Lsn represent, respectively, the mass of the shell, the radius 
and velocity of the shell, the interstellar pressure and its density and the energy release rate 
from an OB association. These equations are, respectively, the equation of motion, the mass 
conservation, the relation of a size R to a velocity v, and the first law of thermal physics. If 
we assume p = (2/3) x E /(4irR 3 /3), these four equations can be solved numerically. Figure 
9 shows a resultant expansion law of a spherical superbubble with Lsn = 3 x 10 37 erg s _1 in 
a uniform interstellar medium of n = 0.3cm -3 . Each curve corresponds to different external 
pressures as p out = 1.7 x lCT 12 erg cm -3 (solid line), p ont = 1 x 10" 12 erg cm -3 (dotted line), 
and p out = (dashed line). Weaver et al.'s (1977) solution, equation (1), agrees with the 
curve of p out = 0. This figure shows that the interstellar pressure plays an important role 
especially in the late phase of the evolution. This is understood as follows: in the late phase 
of the superbubble the difference between the internal pressure and the outer one is small 
and this small difference drives the shell further. Therefore, the critical luminosity would 
be underestimated if we use a solution without taking the outer pressure into account. 
The shell of a superbubble continues to expand as long as the energy ejection continues, 
while a supernova remnant stops its expansion after p ont = p. Thus, exactly speaking, 
the superbubble is never confined as long as the OB association is alive. However, a slow 
expansion driven by a small pressure difference as p — p ont is considered as a signature of 
the confinement. 

In this figure, we also plot an equivalent radius, which is defined using the volume 
occupied with a hot matter (Vhot) as 
("SV \ ^ 3 

-Rcquiv = ( ^ ) j (^9) 

for Models A and C. The equivalent radius for Model A in which the bubble is almost 
confined in the galactic gaseous disk shows a similar expansion law to that obtained by 
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a thin-shell model (a solid curve). In contrast, that of Model C indicates a completely 
different expansion law such that after t ^ 10 Myr the equivalent radius increases rapidly 
and in t ^ 35Myr -R oqu i v surpasses the model with p out = 0. These differences seem to come 
from the distribution of magnetic field strength. In Model A the total pressure (thermal plus 
magnetic one) is almost constant as the bubble expands, because the magnetic pressure is 
dominant over the thermal one and magnetic fields are uniform. While, in Model C the total 
pressure drops according to the density distribution p(z). This figure indicates that when 
the equivalent radius is well fitted by this thin shell model the bubble is nearly confined in 
the disk even if the gas disk has a finite scale-height. In contrast, if hot gas is ejected from 
the galactic disk, the equivalent radius shows a more rapid expansion than that derived by 
this thin-shell model. Values of equivalent radii at the epochs when numerical runs end are 
shown in table 1. 

4.2 Observability 

A shear motion in the galactic rotation and rotation itself may play a role in the evolution of 
a superbubble ( Tenor io-Tagle & Palous 1987; Palous et al. 1990; Silich 1993). Galactic shear 
seems to deform the shape and the Coriolis force makes the shell rotate. The characteristic 
time-scales of the rotation, tr, and the shear, r$ are estimated respectively as 

r R ~ ~ 40Myr(fi /26km s" 1 kpc -1 ) -1 (20) 

and 

r s ~ (IdVl/dR)- 1 ~ 320Myr(Z/lkpc)- 1 (fi /26km g- 1 kpc _1 ) _1 (i2o/8.5kpc), (21) 

where Qq, Rq, and I are the angular speed of galactic rotation, distance from the galactic 
center, and a typical size of a superbubble, respectively. Since active SN-explosion phase 
continues for ~ 50 Myr for an OB association (McCray & Kafatos 1987), in the late phase 
of tr ^ t ^ 50Myr, the effect of the Coriolis force seems to appear as a deformation force 
of the shell. The acu-dynamo mechanism driven by superbubble was studied recently by 
Ferriere (1992). 

The galactic rotation has a little effect on the evolution of a superbubble. Thus, if shells 
or holes observed in external galaxies are elongated after their inclinations are corrected, 
their direction seems to indicate that of the magnetic fields. There have been listed 141 HI 
holes in M31 by Brinks & Bajaja (1986). These holes are observed, more or less, as elliptical. 
Particularly, the holes found near the major axis of M31 are important to determine the 
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physical shape of the holes. Many of these HI holes have their major axes perpendicular to the 
galaxy's major axis. Since this is not explained by projection due to the inclination of M31 
(i = 77 deg), these seem to have physically a shape elongated along the azimuthal direction of 
the galaxy. This is not inconsistent with observations indicating that a global pattern of the 
magnetic field is ring-like in M31, which is measured by radio linear polarization observations 
(for a review, see Sofue, Fujimoto & Wielebinski 1986), in other words, magnetic field lines 
run in the azimuth direction in M31. 

4.3 Porosity 

If hot gas contained in superbubbles occupies a large volume of the galactic disk, a picture 
of the interstellar medium should be changed (McKee & Ostriker 1977). The fraction of 
areas covered by superbubbles younger than r active is estimated with a quantity called as 
two-dimensional porosity which is defined as 

P ^active 

Q(t < r activc ) = r OB / S(t)dt, (22) 
jo 

where tob is the formation rate of OB associations per unit area and S(t) represents the 
area which covered by a hot cavity on the mid-plane of the disk z = 0. This is identical with 
a two-dimensional porosity parameter calculated by Heiles (1990). He estimated galaxy- 
wide average of two-dimensional porosity Q2D — 0.30. r act i ve should be chosen equal to the 
oldest age of a superbubble which contains a hot gas inside. If we assume r act i ve = 20 Myr 
and integrate S(t) for Models A and C (S = irx c y c ), these two models give respectively 
2.08 x 10 6 pc 2 Myr and 1.57 x 10 6 pc 2 Myr. We adopt the estimation of tob from a galactic 
type II supernova rate of rn ~ O.Olyr™ 1 , that is, we assume that all type II SNe occur in OB 
associations, number of type II SNe in an OB association is constant irrespective of richness 
of association as iVsN ~ 100 and OB associations are uniformly distributed in the galactic 
disk with radius i? ga i ~ 10 kpc. This gives an estimation of OB association formation rate 
as 

r OB = N r " 2 , (23) 

iV S N K -f\ a l 

r OB . 3, x lO-Vr" Myr - (^) (^) " (^) 2 . ,24, 

This indicates the two-dimensional porosity to be equal to Q(t < 20Myr) ~ 0.5 — 0.6. Thus, 
rather large fraction of the galactic disk, 1 — exp (-Q) ~ 40% — 45%, is covered by young 
(t < 20 Myr) superbubbles. 
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